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Равновесная плоская трещина с угловыми точками контура в упругом слое 


Б. В. Соболь, Е. В. Рашидова, Е. В. Борисова, С. Б. Петренкова 
(Донской государственный технический университет) 


Рассмотрена трёхмерная статическая задача теории упругости о равновесии упругого слоя, ослабленного 
плоской прямоугольной трещиной. Трещина расположена в срединной плоскости слоя, поддерживается в 
раскрытом состоянии под действием нормальной нагрузки, приложенной к её берегам. Грани слоя находятся в 
условиях гладкого контакта с двумя жёсткими основаниями. Применением двумерного интегрального пре- 
образования Фурье к уравнениям равновесия задача сведена к решению известного сингулярного интег- 
ро-дифференциального уравнения относительно функции раскрытия трещины. Решение уравнения строится 
прямым вариационным методом. В окрестности угловых точек контура решение достраивается численно, с 
учётом заранее выделенной особенности. Получены значения коэффициента интенсивности нормальных 
напряжений в окрестности контура трещины. Установлены особенности поведения решения в окрестности 
прямолинейных участков и угловых точек контура. 
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Введение. Как было установлено в многочисленных работах отечественных и зарубежных ис- 
следователей (например, [1]), в точках гладкости контура трещины функция раскрытия трещины 
должна иметь асимптотику 
у(х, у) =М(х,у) т” (х, у), (1) 
где п(х, у) — расстояние по нормали от внутренней точки области до контура трещины, Мх, у) — 
некоторая регулярная функция. 

В случаях, когда контур трещины имеет особые точки: во-первых, нарушается состояние 
плоской деформации в окрестности этих точек, во-вторых, меняется показатель особенности поля 
напряжений и перемещений в окрестности контура вблизи таких точек. Определению показателей 
при особенности в зависимости от угла излома контура посвящён ряд работ отечественных и за- 
рубежных исследователей [2]. 

Если в плоскости расположения трещины ввести полярную систему координат с полюсом в 
особой точке контура, то функция раскрытия трещины в окрестности этой точки асимптотически 
может быть, описана соотношением 


Убиу)=Ца) [8], © 


где р, а — полярные координаты точки, [(а) — некоторая регулярная функция. 

Рассмотрим наиболее распространённый частный случай, когда контур трещины в точке 
излома образует прямой угол. В этом случае установлен следующий показатель особенности в 
соотношении (2): = = 0,816. 

Постановка задачи. Продемонстрируем методику решения задач теории упругости для трещин с 
особыми точками контура на примере задачи о прямоугольной трещине в слое. 

Трещина расположена в срединной плоскости слоя, поддерживается в раскрытом состоянии 
под действием нормальной нагрузки, приложенной к её берегам. Грани слоя находятся в условиях 
гладкого контакта с двумя жёсткими основаниями. Применением двумерного интегрального пре- 
образования Фурье к уравнениям равновесия в перемещениях задача о продольной равновесной 


60 


Вестник ДУ, 2012. № 5 (66) 








трещине в слое сводится к решению двумерного сингулярного интегро-дифференциального 
уравнения [3]: 


АУ (Е, ие + [у(Е.п)5(х, у, п) ай = При у), (3) 





здесь (хуеО, О — область расположения трещины в плане, А = Ух -&) +(у -п), 
2. 2 

9=Е / [2 (1 - У) | ‚А= р. -- р Б(х, у) — нагрузка, приложенная к берегам трещины, Е — модуль 
Хх у 

Юнга, их — коэффициент Пуассона, $(х,у,Е,п) — регулярная часть ядра уравнения, характе- 


ризующая влияние геометрических параметров задачи: размеров трещин в плане, относительных 
расстояний между ними и (или) до границ рассматриваемых тел. Для задачи о трещине, лежащей в 
срединной плоскости слоя, имеем [3]: 


бу п) = (В В), Е (@) = [ш[н(и)- и] (аи, 


Н(и) =и ($12и +2и) (спи Ен 1 (х) — функция Бесселя 1-го рода. 
Решение уравнения (3) будем строить прямым вариационным методом, описанным в [4]. В 
рассматриваемом случае область интегрирования описывается неравенствами: © :|х|<а, |у|<В 


(рис. 1). 





Рис. 1. Схема трещины 


Метод решения. Выберем, систему координатных функций таким образом, чтобы установленная 
выше асимптотика (2) с соответствующим показателем явно учитывалась. Ограничимся рассмот- 
рением, как наиболее характерного частного случая, равномерно распределённой нагрузки 
р(х,у) на берегах трещины, где р(х,у) = р = сопзЁ . Решение задачи при этом строится в виде: 


чех] 





М 
у(х, у) = ва" (х, уу Ат соз ИХ соз ПХ (4) 
т=0 п=0 а р 
Коэффициенты А» в разложении и определяются из условия минимума функционала: 
У) = Пу бк,у) ДУ &1)К, (Е -х,п- узаванан, 717 [уху рава, (5) 





где К, (а, В) = Пл ИМ +2 }созаи со5Видиам . 
00 


Следует отметить, что в аналогичном виде вариационный метод был использован при ис- 
следовании задачи о равновесии упругого пространства ослабленного трещиной [5]. В качестве 
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координатных функций в этой работе были использованы степенные функции с множителем, явно 
содержащим требуемую асимптотику в окрестности контура трещины. 

Минимизация функционала (5) в рассматриваемом случае приводит к следующей системе 
линейных алгебраических уравнений для определения коэффициентов А): 


м м 
У А.О т =27 (пт,пп), (6) 
1=0 7=0 
где т =0,1,2,...М; п=0,1,2,....М. 
7, (а)./7 
В данной задаче гв) ОВ. При вычислении коэффициентов системы (6) ис- 


пользовано значение интеграла 


2 
Г 41? (х,у) 0$ 4 ахау = ПавЕ (а,В), 
К а Ь 4 





полученного с помощью (3.752) из [6]. 
Вычисление элементов четырёхмерной матрицы О., производится с использованием 


уп 

представлений (5.4) из [4]. 

Сравнительный анализ и сопоставление численных результатов. Конечной целью иссле- 
дования задачи теории упругости для тела с трещиной, как известно, является определение кри- 
тических нагрузок, при которых трещина становится неустойчивой и начинает расти. Поскольку в 
качестве основного критерия разрушения в данной работе принят силовой критерий [7], перейдём 
к определению коэффициента интенсивности нормальных напряжений в окрестности контура 
трещины в каждом из рассмотренных случаев. Особое внимание в дальнейшем также будет уделяться 
сравнительному анализу применяемых здесь методов и сопоставлению численных результатов. 

Контроль точности решения задачи вариационным методом приводился сопоставлением 
результатов при различном количестве используемых координатных функций (Ми М. В частности, 
вычислялась амплитуда раскрытия берегов трещины в центре при различных соотношениях гео- 
метрических параметров задачи. Во всех рассмотренных ниже случаях расхождение результатов 
счёта при М =М =4 и М =М =6 не превысило 3,5 %. 

Кроме того, о достоверности получаемых результатов могут свидетельствовать и следую- 
щие качественные рассуждения. При вытягивании трещины вдоль одной из осей симметрии, ам- 
плитуда её раскрытия в центре, а также численные значения коэффициента интенсивности нор- 
мальных напряжений в средней части большей из сторон прямоугольного контура должны быть 
близки к известным соответствующим решениям задач о полосовой трещине и эллиптической 
трещине с такими же полуосями. Причём, поскольку упомянутая полосовая трещина ширины 26 
содержит внутри себя прямоугольную область, а эллиптическая, в свою очередь, содержится в 
последней, решение задачи о прямоугольной трещине, по известному принципу сравнения [8], 
должно быть заключено между соответствующими решениями этих двух задач. В таблице 1 пред- 
ставлены приведённые значения амплитуды вертикальных перемещений у (0,0)6 / (ра) в центре 
трещины для полосовой [3], эллиптической [9] и прямоугольной трещины. Вычисления проведены 
для предельного случая Л -> о . Для сопоставления, здесь же приведён соответствующий результат 
работы [5]. 

В случаях эллиптической и прямоугольной трещин при вычислениях принято следующее 
соотношение размеров Б /а=0,5. Вычисления показывают, что и при конечной толщине слоя 


установленные закономерности не нарушаются. 
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Таблица 1 
Приведённые значения амплитуды вертикальных перемещений в центре полосовой, 
прямоугольной и эллиптической трещины 


(4) [5] [9] [3] 
0,882 0,872 0,826 1,000 


Перейдём к анализу численных значений коэффициента интенсивности нормальных на- 
пряжений К, при обходе по контуру прямоугольной трещины. В этой задаче более удобно рас- 


сматривать не относительные значения параметра К; (как в задачах об эллиптических трещинах, 


где точное аналитическое решение для предельного случая Л -> х существует), а абсолютные. 
Вдоль каждой из сторон прямоугольной области значения коэффициента интенсивности нор- 
мальных напряжений будем вычислять по формулам: 

тп 


2 
К» (х) = ра 15! (х,0) ‚> - 


2 мм 
Кь(У) = ра =1(0,У) > > Асов. (8) 
т=0 п=0 


Графики изменения коэффициента интенсивности нормальных напряжений вдоль большей и 
меньшей сторон прямоугольника представлены на рис. 2 (а, 6). Вычисления проведены для ряда 
значений параметра Л =А/а при Б /а =0,5. Отметим, что максимальные амплитудные значения 


М 
УА т, ый (7) 


[М 


З 
И 


К, достигает на большей стороне. Опираясь на силовой критерий разрушения, можно сделать 


вывод, что при достижении приложенной нагрузки своего критического значения, рост прямо- 
угольной трещины начнётся с середин больших сторон. На рис. 3 представлены соответствующие 
графики, иллюстрирующие решение задачи для случая трещины квадратной формы (Ь = а). 

С приближением к угловым точкам контура трещины коэффициент интенсивности нор- 
мальных напряжений стремится к нулю. Это обусловлено особым видом асимптотики (2) функции 
раскрытия трещины в окрестности этих точек, с соответствующим показателем. Перейдём к по- 
строению решения задачи у(х, у) в окрестности угловых точек. 

Вдоль лучей, выходящих из угловой точки контура под различными углами а, вычисляются 
значения функции: [ (а) =у(х,у)/(р/а). 

Фрагменты результатов таких вычислений представлены в таблице 2 (при Ма = 0,5) и 
таблице 3 (при Ва = 1). 

В каждом из случаев, представленных в таблицах, нетрудно обнаружить интервал изме- 
нения отношения р/ ак Ме ‚ в котором функция с достаточно высокой точностью сохраняет 


стабильность вдоль каждого луча и изменяется лишь при переходе с одного луча на другой. Упо- 
мянутая закономерность ранее была установлена для предельного случая рассматриваемой здесь 
задачи, соответствующего Л-›со [5], [10]. 

Расположение и размеры этой зоны зависят от а, отношения сторон области трещины Б/аи 
относительной толщины слоя Л. Таким образом, определив для каждого значения величину [(а) на 
указанном отрезке [м, >], представляется возможным достроить искомую функцию у(х, у) по 
асимптотической формуле (2) непосредственно до вершины угла трещины. 

Графики функции (а) при различных относительных толщинах слоя представлены на ри- 
сунке 4 (слева — Б/а = 0,5, справа — Б/а = 1). 

Следует отметить симметричность функции у(х, у) в окрестности угловой точки контура 
относительно биссектрисы. Между тем, сама прямоугольная область расположения трещины такой 
симметрией не обладает. 
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0 0,25 0,5 0,75 х/а 





0 0,33 0,66 х/а 


Рис. 2. Изменение коэффициента интенсивности нормальных напряжений вдоль 
большей (а) и меньшей (6) сторон прямоугольной трещины 
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Рис. 3. Изменение коэффициента интенсивности нормальных напряжений для случая трещины квадратной формы (В = а) 


Таблица 2 
Приведённые значения функции раскрытия вдоль лучей, выходящих из угловой точки 
прямоугольной трещины (Б/а = 0,5) 
1/2 


Таблица 3 
Приведённые значения функции раскрытия вдоль лучей, выходящих из угловой точки 
квадратной трещины (Б/а = 1) 





уе 
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М 


а) 6) 


40 
Рис. 4. Изменение множителя — (а) при особенности функции раскрытия трещины в окрестности 
ра 





угловой точки для прямоугольной (а) и квадратной (6) формы трещины 


Выводы 

1. В результате проведения численного эксперимента установлено, что решение задачи о 
прямоугольной трещине в упругом слое удовлетворяет принципу сравнения. 

2. Коэффициент интенсивности напряжений достигает своего максимального значения при 
обходе по контуру в окрестности центра длинной стороны прямоугольной области. 

3. Установлена симметричность функции раскрытия трещины в окрестности угловой точки 
контура, независимо от соотношения размеров прямоугольной области. 
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ЕОЧТЫВВТОМ РЕАМАК СКАСК МТТН СОМТОЧК АМСУГАК РОТМТ$ 
ТМ ЕСАЗТТС ГАУЕВ 


В. М. боБо/, Е. \. ВКазМаома, Е. \. Вом5ома, $. В. РегепКо\а 
(Роп Зае Тесптка! Упмег$Ку) 


Тре гее-айтепхюпа! “айс ргоМет ое ейзйсйу {теогу оп {пе е/а5Нс !ауег Ба!апсе игеакепеа Бу а ПаЁ гесапдишаг 
сгаск #5 сопз!егеа. Тре сгаск 15 осайеа т пе тебап 1ауег риапе; # /5 зирропеа ореп ипаег пе погта! са тд 
арр/еа Ю #5 еддез. Тре !ауег ед дез аге ипаЕег {те зтооЁ? сопёасЕ илёй о поГа Базез. ТНе рго ет /5 сопуееа ю 
{Пе зо/ивоп оЁ {пе Кпоит ташаг и®едго-@егепва! едиайоп сопсегптд пе сгаск орептд ипсвоп гоидй пе 
арр/сайоп о! вмо-@тепзюопа! итедга! Роипег капзГогтайоп 0 те едийЬпит едиайоп. Тйе едиайоп зо/иНоп 15 
сопзтисеа Бу {Пе СгесЕ уапавопа! теёпод. тп пеоПБогпооа о! {пе сотоиг апдшаг ротЁз, йе о/ивоп сотр/еЁе5 
{пе сопзбгисйоп питепсайу, ий гедага Гог {пе еапу @5Нпошевеа этдашат у Уаиез оЁ {пе @гесЕ ${те5 п{еп5Йу 
гасюг т педИБогптоо4 ог пе сгасК репрВеГу аге оббеатеа. Тре 5о/ийоп Берауюга! ЕеаЁигез т пеовБогрооа о! {те 
сотюиг гесШтеаг $йез апа апдиш!аг рот аге еза вед, 

Кеуигогаб: сгаск, сотоир 5тез$ И{еп5йу асю; (вуер тдашаг Иедга! едиайоп, иапайопа/ те под, Гогсе басиге 
стШепоп, за!епё ротё пеНБогптоод, !ауег геаН\е Е/СКпеЕз5. 
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